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Abstract 

This paper investigates correlation properties of fluctuations in fatigue crack growth of 
polycrystalline materials, such as ductile alloys, that are commonly encountered in structures 
and machinery components of complex electromechanical systems. The model of crack damage 
measure indicates that the fluctuations of fatigue crack growth are characterized by strong 
correlation patterns within short time scales and are uncorrelated for larger time scales. The 
two correlation regimes suggest that the 7075-T6 aluminum alloy, analyzed in this paper, is 
characterized by a micro-structure which is responsible for an intermittent correlated dynamics of 
fatigue crack growth within a certain scale. The constitutive equations of the damage measure 
are built upon the physics of fracture mechanics and are substantiated by Karhunen-Loeve 
decomposition of fatigue test data. Statistical orthogonality of the estimated damage measure 
and the resulting estimation error is demonstrated in a Hilbert space setting. 
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I. INTRODUCTION 



The fracture of solids and the growth of cracks is a typical instability phenomena which 
are known to be strongly nonlinear. Herein we apply to fracture mechanics some of the recent 
methods developed in statistical physics. In particular, we use the notion of fractal statistics 
to describe the correlation of the fluctuations around fatigue crack growth in polycrystalline 
materials, such as ductile alloys. In this paper, we have investigated the fatigue fracture 
properties of 7075-T6 aluminum alloy. 

The importance of this investigation is that, in both the physics and engineering literature, 
the fluctuations around fatigue crack growth in a typical material have always been assumed 
to be random or uncorrelated noise. Consequently, the associated models include uncorre- 
lated random processes. For example, in agrement with the existent theory of micro-level 
fatigue cracking, Bogdonoff and Kozin proposed a Poisson-like uncorrelated-increment 
jump model of fatigue crack phenomena. An alternative approach to stochastic modeling of 
fatigue crack damage is to randomize the coefficients of an existing deterministic model to 
represent material inhomogeneity j^J]. A third approach has been to adopt a deterministic 
model of fatigue crack growth in addition to a random process, see for example 

The fatigue crack growth process can also be modeled by nonlinear stochastic differen- 
tial equations using It 6 statistics [fl that again presuppose randomness of the fluctuations. 
Specifically, the Kolmogorov forward and backward diffusion equations, which require solu- 
tions of nonlinear partial differential equations, have been proposed to generate the statis- 
tical information required for risk analysis of mechanical structures [?], ls[ • These nonlinear 
partial differential equations have only been solved numerically and the numerical proce- 
dures are computationally intensive as they rely on fine-mesh models using finite-element 
or combined finite-difference and finite-element methods [3]. Casciati et al. have ana- 
lytically approximated the solution of the ltd equations by Hermite moments to generate a 
probability distribution function of the crack length. 

Several studies have determined that the stochastic fluctuations observed in innumerable 
natural phenomena are not simply random, that is, uncorrelated noise, but present correla- 
tion patterns that reveal complex and alternative dynamics and/or material microstructures. 
Thus, the purpose of the present research is to determine whether uncorrelated stochastic 
models such as those previously discussed in the literature are realistic in describing the 
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fluctuations around fatigue crack growth in polycrystalline materials, or whether such fluc- 
tuations present patterns that would reveal complex material micostructure requiring alter- 
native correlated stochastic models. Two main classes of correlation patterns are commonly 
observed in natural time series and these are denoted as short and long-time correlations. 
Short-time correlations are characterized by phenomena that rapidly lose memory of past 
or distant events. This happens, for example, when the autocorrelation function of the time 
series decays exponentially in the time separation between two elements. By contrast, long- 
time correlations are characterized by autocorrelation functions that decay more slowly than 
(negative) exponentials; one example is the inverse power-law decay. 

A simple model, which has been extensively used in the interpretation of stochastic 
fluctuations in a time series with % = 1,2, ... ,N, is based on the evaluation of the 
mean-square displacement of the diffusion-like processes generated by trajectories X n (t) 
defined as 



If {C,i} is a white random sequence, the diffusion process is a well-known Brownian motion. 
The central limit theorem applied to the diffusion distribution generated by trajectories 
X n (t) yields a probability density that converges to a Gaussian function whose mean-square 
displacement converges asymptotically to 



with a = 1. In general, it is possible to have anomalous behavior yielding enhanced diffusion 
(a > 1) that has been known for twenty years to arise in dynamically chaotic systems 11 1, 
or sublinear diffusive growth (a < 1) that is familiar from disordered fractal materials 12 1. 

Anomalous diffusion reveals persistent (for an enhanced diffusive growth) or antipersistent 
(for a sublinear diffusive growth) correlation patterns in the dynamics of a random walk. A 
persistent random walk is characterized by a probability of stepping in the direction of the 
previous step that is greater than that of reversing directions. An antipersistent random 
walk is characterized by a probability of stepping in the direction of the previous step that is 
less than that of reversing directions. Sometimes a momentarily initial enhanced or sublinear 
diffusive growth, lasting up to a certain time-scale, is generated by the statistical transition 
to the asymptotic regime of the diffusion process. For example, a simple discrete random 
walk is described by a binomial distribution that only asymptotically converges to a Gaussian 
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while initially presenting an enhanced diffusive growth Q|. Thus, a real autocorrelated time 
series will lose its correlation patterns if the temporal order of the sequence is randomized. 

There are a number of different theoretical approaches that explain the anomalous dif- 
fusion depicted in (J2J). One such explanatory model is that of an infinitely long correlated 
random walk in which a = 2H, where H is the Hurst exponent in the interval < H < 1 
with the case H = 0.5 corresponding to a simple random walk. This model has been used 
extensively in the interpretation of fluctuations in time series in the physical and life sciences 
\h\ and is called fractional Gaussian noise Q| . Another kind of anomalous diffusion has to 
do with taking steps that are uncorrelated in time, but on a random or fractal, not a regular 
lattice. In the second model, an anomalous diffusion occurs because geometrical obstacles 
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exist on all length scales and such obstacles inhibit transport. Havlin and Ben-Avraham 
point out that the anomalous exponent a is related to the fractal dimension of the random 
walk path on the lattice. There is a third possible explanation of the anomaly in @ called 
a Levy walk that was first used to understand turbulent diffusion Q| and yields a ~ 3, 
which is consistent with Richardson's law of enhanced diffusion ^1 . 



Physical examples of anomalous diffusion processes are earthquakes 



tur 



)U 
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L rainfall 
2111 , solar flares |2J 



ent fluid flow |20j, relaxation of stress in viscoelastic materials [1J, 
24j , and other processes with slip-stick dynamics. Recently, a multi-scaling comparative 
analysis to distinguish Levy walk intermittent noise from fractal Gaussian intermittent noise 
was suggested by Scafetta and West j^ . 

Finally, a physical system might be characterized by different values of the scaling ex- 
ponent a at different scales j^J. Usually, this means that one system is characterized by a 
non-self affine structure. The scale at which the transition from a scaling regime to another 
occurs indicates the scale at which the structure changes. In this work we determine that 
the fluctuations around the ballistic growth of fatigue cracks in ductile alloys present such 
a scale transition from a strongly correlated regime at short-time scales to a random regime 
at longer time-scales. Properties, such as grain size distribution, degree of heterogeneity, 
the existence of microscopic defects, inclusions, twin boundaries and dislocations, of poly- 
crystalline materials may contribute to the micro-mechanisms of fatigue fracture revealed 
by the present analysis. 

This paper is organized into six sections, including the present one, and an Appendix. 
Section |H] provides the underlying phenomenology of the stochastic damage measure. Sec- 
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tion IIIII presents Karhunen-Loeve (KL) decomposition of fatigue test data to formulate an 
estimate of the stochastic measure, which is statistically orthogonal to the estimation error. 
Section HVI focuses on identification of the model parameters and their probability distribu- 
tions. Section |3 presents the results of model prediction by Monte Carlo simulation. The 
paper is summarized and concluded in Section lVTl with recommendations for future research. 



II. MEASURE OF FATIGUE CRACK DAMAGE 



Traditionally fatigue crack growth models have been formulated by fitting estimated mean 
values of fatigue crack length a t , generated from ensemble averages of experimental data, 



as functions of time in units of cycles 



221, & Ray and Patankar 



have formulated the 



state-space modeling concept of crack growth based on fracture-mechanistic principles of the 
crack-closure concept j^. The state-space model has been validated by fatigue test data 
for variable-amplitude cyclic loading, see for example Refs. j^sl. I31I l^. 

The three panels in Figure 1 show test data of cumulative fatigue crack growth in the 



33| . It is important to note that 



7075-T6 aluminum alloy under different cyclic loading 
the crack growth curves do not increase smoothly, but they exhibit fluctuations around an 
ideal smooth curve of crack growth representing ballistic growth. In this context, a major 
objective of the paper is to investigate the autocorrelation properties of these fluctuations 
with the smooth curve removed. In the following we briefly review the theory and the 
standard phenomenological equations that describe the fatigue crack growth. 

In linear fracture mechanics, it is assumed that the stressed material remains elastic 
and undamaged everywhere, except in a small domain in the vicinity of the crack tip. 
However, this view is not confirmed by experimental evidence and the process of fatigue 
damage accumulation could occur throughout the stressed volume. Paris and Erdogan 27 1 
originally developed a phenomenological model of crack growth rate, which depends on the 
stress history and is thus represented by a continuum rate equation having the hereditary 
structure. This model has been subsequently modified by many researchers (see, for example, 



citations in Refs. 



34, 
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in the following form. 
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t = Q>t — at- st 



h (aK } 



eff 



St, 



(3) 



with h(0) = and a to > for t > to, where a t is the estimated mean of the crack length at 
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the time t during a stress cycle and St is the time duration of the stress cycle; and AK\^ 
is the stress intensity factor range at time t, which is given by the experimentally validated 
empirical model. 

AK? f = AS^tt atst F(a t _ st ), (4) 

where ASt is the range (i.e., the difference between maximum and minimum values) of the 
stress cycle at time t, which is directly related to the applied load. Experimental observations 
suggest that both duration and shape of a stress cycle are not relevant for crack growth in 
ductile alloys at room temperature. A stress cycle is only characterized by the minimum 
stress S min and the maximum stress S max , respectively, and is denoted as the ordered pair 
^gmin^gmaxy rpj^ em pj r j ca j re l a tion .F(«) in Eq. (j3J) represents the geometry of the crack 
tip; for center-cracked specimens of half- width w with < at < w at all t > to, the structure 
of F(») has been experimentally determined as l^ : 
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F(a t st) = ^sec l^— a t ^tj- (5) 

The function h(») in Eq. Q is a non-negative Lebesgue-measurable function that is 
dependent on the material and geometry of the stressed component. It has been shown in 
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that, for center-cracked specimens of ductile alloys, 



the fracture mechanics literature 
the function h(») obeys the power law: 

fc(AJS // ) = (A*?T> ( 6 ) 
where the exponent parameter m is dependent on the material of the stressed component; 
for ductile alloys, m is in the range of 2.5 to 5.0 j^]. 

Equations Q, (JH), © an d © ar e now combined to formulate a mean- value model of 



fatigue crack growth for center-cracked specimens of ductile alloy materials: 

r / 1 m 

5a t oc AS t \\h-st sec ( ^ atst) St. 



(7) 



with a tn > and t > tn. 

a 

Following Sobczyk and Spencer [2( and the pertinent references cited therein, we random- 
ize the deterministic mean- value model, Eq. (jJJ), to obtain a stochastic model for the rate 
of crack growth. The stochastic model of continuous crack length is built upon the model 



structure proposed by Ray 
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dct(0 = n((,t) 



ASt 



and is given by: 

i 

dt 



ct(C) 



^cos(f c t (0 



n((,t) 



(ASt ygg 

l-m(f ( 



dt. 



(8) 
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where the random sample ( signifies a specimen or a machine component on which a fatigue 
test is conducted; the dimensionless stochastic crack length Ct(C) is normalized with respect 
to the half width w, that is, the mean value c t = a t /w. Equation JHJ) is a continuous 
stochastic version of Eq. (J7J), where the differential of the stochastic crack length dc t (() is a 
function of the crack length c t (() at time t and the normalized stress AS t = ASf/S y , where 
S v is the yield stress of the material. The condition < Ct < q < is imposed to ensure 
non- negativity of the crack length increment almost surely, i.e., dc t (() > for almost all 
samples (. The stochastic process of crack growth is largely dependent on the second-order 
random process Q((,t) and the exponent parameter m in Eq.JHJ). 

To investigate the stochastic properties of the fatigue crack growth process, we separate 
fi(C, t) into two parts as: 

n(C,t) = n (C)[i + ni(C,*)], (9) 

where the time-independent component Qq (() represents uncertainties in manufacturing, 
for example in machining, and makes a major contribution to the ballistic component of the 
crack growth; the time-dependent component Q± {(,t) represents uncertainties in the mate- 
rial microstructure and crack length measurements that may vary with crack propagation in 
a sample (. This latter component is primarily responsible for the small fluctuations around 
the ballistic component of crack growth whose autocorrelation properties we study. 

We postulate that Qo and Q\ in Eq. Q are statistically independent of one another for 
all t > to, where to is the initial time. The rationale for this independence assumption is 
that inhomogeneity of the material microstructure and measurement noise, associated with 
each test specimen and represented by fli (C,£), are unaffected by the uncertainty Q (() 
due, for example, to machining operations. Without loss of generality, we assume that the 
fluctuations in time have a zero mean value, i.e., (Qi {(,t)) = for all t > t . Furthermore, 
non- negativity of the crack growth rate dct (C) in Eq. (jHj) is assured in the almost sure (a.s.) 
sense by imposing the constraint Qq (() > with probability 1 (wp 1). 

For notational brevity, let us suppress the term ( in random processes like c t (() and 
Q (C, t). A combination of Eqs. (jHJ) and © and few simple algebraic steps yield the following 
equation for each sample point £ : 

dc t = {AS t ) m tt [1 + fti (t)} dt w.p.l. (10) 
Pointwise integration of Eq. (fTU|) yields the solution of fatigue damage increment from the 
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initial time t Q to the current time t as: 



tf, (t, t ) = J (AS t ,) m n [1 + fli (f ) } dt' w.p.l 

to 



(11) 



An explicit expression of the stochastic diffusion process ip (t, to) is obtained by integrating 
the left side of Eq. ()1L)|) and is given by 



iP(t,t ) = 



l-m/2 



l-m/2 



m/2 



m 



3-m/2 



3-m/2 
C <0 



3 -m/2 



(12) 



where ip (t, t Q ) represents a dimensionless non-negative measure of fatigue crack damage 
increment from the initial instant to to the current instant t as a function of the normalized 
crack length. The constant parameter m in (J 12)) is in the range of 2.5 to 5 for ductile alloys 
and metallic materials ensuring that (1 — m/2) < and (3 — m/2) > 0. The diffusion process 
ip (t, t ) is almost surely continuous because it is a continuous function of the crack length 
process c t wp 1. Both c t and x/i (t,t ) are measurable functions although their (probability) 
measure spaces are different. In essence, the probability of ip(t,to), conditioned on the 
initial crack length q , leads to a stochastic measure of fatigue crack damage increment at 
the instant t starting from the initial instant to- 

For a constant stress range AS", we carry out the time integration in Eq. fTTf to obtain 



V (t, to) = (AS) m [fi (t - to) + 8(t, to)] 
where the second term on the right side is the time integral 
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e(Mo) = n [^i (t')dt'. 



(13) 



(14) 



Thus, the stochastic diffusion process ip (t, to) according to the model (|T3*|) is given as the 
sum of a random component, linear in time, plus a time-fluctuating component proportional 
to the diffusion process (t, t ). 

The objective is to validate the model in Eq. (|11|) by decomposing the damage increment 
measure i/j (t, t ) into two parts that are mutually statistically independent and, at the same 
time, equivalent to the two components of the right side of Eq. (|13j). That is, we would like 
to obtain an estimate ip (t, to) of the stochastic damage increment measure ip (t, to) and of 
the fluctuations i/j (t, to) around i/j (t, to) from the initial instant to to the current instant t 
such that: 

V>(Mo) = ^(Mo) + <MMo) , (15) 
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where ip (t, to) is statistically equivalent to AS m flo(t — t ), and z/> (t, t ) is statistically equiv- 
alent to (AS) m Q{t,t ) of Eq. ([T31) . 

To test the validity of the above postulate that the two components of the multiplicative 
random process Qq (£) and ((,t) in Eq. Q are statistically independent, we require that 
the zero-mean estimation error ip (t, to) be statistically orthogonal to the estimate of the 
increment measure ift (t,to) in the Hilbert space L2 (P) defined by the probability measure 
P. As such if; (t, to) is the best linear estimate of the stochastic diffusion process. Based 
on mean-square continuity of the damage measure ip (t, to), the next section elaborates on 
the model structure laid out in Eq. (J15|) . To this end, we ana lyze experimental data sets 
of random fatigue via Karhunen-Loeve (KL) decomposition 38, 39, 4^ that guarantees the 
above statistical orthogonality among the components of the decomposition. In Section IV 
we also use these experimental data sets to identify the model parameters. 

III. KARHUNEN-LOEVE DECOMPOSITION OF EXPERIMENTAL DATA 

In this section we analyze fatigue test data via KL-decomposition j^j] to justify the model 
structure postulated in Eqs. (JTTJ) and (fT2|) . We use the experimental data of random fatigue 
crack growth in the 7075-T6 aluminum alloy j^] and conduct the tests under different 
constant load amplitudes at ambient temperature. For all experiments the half-width is 
w = 50.8mm, the initial crack length is a to = 9mm, and, therefore, the initial dimensionless 
crack length is c to = a to /w = 0.18 with probability 1. The Ghonem data sets were generated 
for 60 center-cracked specimens each at three different constant load amplitudes: (i) Set #1 
with peak nominal stress of 70.65 MPa (10.25 ksi) and stress ratio R = S mm /S max = 0.6 
for 54,000 cycles, the effective stress range A5" e = 15.84 MPa; (ii) Set #2 with peak nominal 
stress of 69.00 MPa (10.00 ksi) and R = 0.5 for 42,350 cycles, and AS e = 17.80 MPa; and 
(iii) Set #3 with peak nominal stress of 47.09 MPa (6.83 ksi), R = 0.4 for 73,500 cycles, 
and A5* e = 13.24 MPa. The three experimental datasets are shown in the three panels 
of Figure 1. 

The KL-decomposition requires the mean and covariance of the stochastic measure of 
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damage increment ip (t, to) which are expressed as: 

m (t,t ) = (ip (t,t )) 

(16) 

C^p (t 1 ,t 2 ;t Q ) = ([^ (tx,t Q ) - ^ (ti,t )][il) (t 2 ,t ) - n^> (h,t )}) 

The covariance function C^(ti, t 2 ; t Q ) in Eq. (fTB^) is continuous at t\ — t 2 — t for all t > t . 
Hence, the process ip ft, tn) is mean-square (ms) continuous based on a standard theorem of 



mean-square calculus [38|, |39( . The mean and covariance are calculate for the 60 available 
center-cracked specimens in each case. 

Since only finitely many data points at n discrete instants are available from experiments, 
an obvious approach to the analysis of the damage estimate is to discretize over the finite 
time horizons [to, t] so that the stochastic process ip (t, to) becomes the n-dimensional random 
vector ip. Consequently, the covariance function C^(tx,t 2 ;t ) in Eq. (fTB|) is reduced to a 
real semipositive-definite (n x n) symmetric matrix C^. Since the experimental data were 
collected at sufficiently close intervals, contains pertinent information of the crack 

damage process. The n (real non-negative) eigenvalues of are ordered as Ai > A2 > 
• • • > A n , with the corresponding eigenvectors, y? 1 , if 2 , • • • , <p n , that form an orthogonal 
basis of 3?" for signal decomposition. The KL-decomposition also ensures that the n random 
coefficients of the basis vectors are statistically orthogonal, that is, they have zero mean 
and are mutually uncorrelated. These random coefficients form a random vector X = 
[x\ ,x 2 , ■ ■ ■ ,x n ] having the covariance matrix Cxx = diag (X±, X 2 , • • ■ , X n ) leading to 
a decomposition of the discretized signal as: 

V>=W+X>^' (17) 
j'=i 

It was observed by Ray (3?! that the statistics of crack length are dominated by the random 
coefficient corresponding to the principal eigenvector (i.e., the eigenvector associated with 
the largest eigenvalue) and that the combined effects of the remaining eigenvectors are small. 
Therefore, the signal ip in Eq. ()17j) is expressed as the sum of a principal part and a (zero- 
mean) residual part that are mutually statistically orthogonal: 

1 

V>= W + + ( 18 ) 

principal part i=2 

residual part 
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Thus, as Eq. (|T5|) requires, the vector tp is expressed as the sum of the principal and residual 
parts with equality in the mean square (ms) as: 



if> = $ + i) (19) 

where the principal part is the damage estimate 

$ = (if;) + xkP 1 , (20) 

the residual part is the estimation error representing the fluctuations around the mean 
damage estimate (j20|) 

i = X>^'> (2i) 

" 2 n 

and the resulting (normalized) mean square error |40|] is: 



s 2 



Trace [Cov ib — $ \ J-L ^ 

J/ - J=2 (22) 



Trace {Cov [if,]} £ A . 



The KL-decomposition of fatigue test data sets reveals that 0.01 < e 2 rms < 0.1 for all three 
data sets. 

The principal eigenvector 1 (t), associated with the largest eigenvalue Ai, closely fits the 
ramp function (t — t ) for each of the three data sets in Figure 1; this is shown in Figure 2 for 
the data set 1. Comparing the terms on the right hand side of the discrete model in Eq. (119(1 
with those of the continuous model in Eq. 1)13(1. it is reasonable to have the random variable 
AS m [Qq — /io] equal (in ms sense) to the random coefficient X\ of the principal eigenvector 
tp l (t). Applying the lemma from the Appendix, a mean-square equivalence between the 
KL-decomposition model in Eq. ()19|) derived from the test data and the postulated model 
in Eq. (fT7)) is established as: 

m)) ™ S AST ft (t- to) (23) 

discrete model (test data) continuous model (constitutive relation) 

xitfit) w AS m [n - /i ] (t - to) • (24) 

discrete model (test data) continuous model (constitutive relation) 

_ n _ mq 

« (AS) m e(t,t ) (25) 

i=2 continuous model (constitutive relation) 

discrete model (test data) 
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Thus, we have $ = (i/j(t)) + xrf 1 (t) w AS m Q (t - t ), and ^ = £ x^' pa AS m 6 (t, i ) 
as assumed in Eq. (|15j) . The two entities on left hand side in Eqs. f!24j) and (|25jl are 
mutually statistically orthogonal by construction. Similarly, in view of Eq. (|15[). the zero- 
mean estimation error ip (t, t Q ) is statistically orthogonal to ip (t, to) in the Hilbert space 
L2 (P) defined by the probability measure P associated with the stochastic process ip (t, to). 
As such ip (t, to) can be viewed as the best linear estimate of ip (t, t ) with the least error 
ip (t, t ) in the mean-square sense. 

It follows from Eqs. (fTHj) to (|25jl that the uncertainties associated with an individual 
sample resulting from the damage measure estimate ip (t, to) dominate the cumulative ef- 
fects of material inhomogeneity and measurement noise in the estimation error ip (t, to) unless 
(t — to) is small. Therefore, from the perspectives of material- health monitoring, risk anal- 
ysis, and remaining life prediction where the inter-maintenance interval (t — to) is expected 
to be large, a reasonably accurate identification of the mean /io and variance of the ran- 
dom parameter Qq is crucial, while the role of the diffusion process (t, to) is relatively less 
significant. This observation is consistent with the statistical analysis of fatigue test data by 
Ditlevsen where the random process described by Eq. (|25|) is treated as the zero-mean 
residual. Ditlevsen [^] also observed largely similar properties by statistical analysis. Never- 
theless, the stochastic properties of fluctuating function 9 (t, to), which we investigate, can 
disclose important information about the material structure of alloys during crack damage. 



IV. DATA ANALYSIS 

In this section we investigate the stochastic equivalence made in Eq. (|23|) between the 
residual component of the signal as obtained by the KL-decomposition and the linear ap- 
proximation. The first step is to evaluate the exponent parameter m by fitting the data of 
the crack growth with Eq. (jHJ). The fit is done by considering the crack increments from all 
60 cases for each of the three experiments. 

By using the empirical values of m it is possible to estimate ip(t,to) via Eq. (|12|). The 
three plots in Figure 2 compare the curve ip(t,to), its principal part according to the KL- 
decomposition and its linear approximation according to the continuous model made in 
Eqs. (|2~3j) plus (j21| for set #1: the figures for the other data sets look qualitatively similar. 
Figure 3 shows the quality of the equivalence made in Eqs. (}2"3*j) plus (j2H) between the 
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discrete model, which makes use of the KL-decomposition, and the continuous model, which 
makes use of a linear approximation. 

Figure 3 shows the fitted data and the results for set #1; the figures for the other sets 
are similar. The parameters for all three sets are listed below. 

• D, (C) AS m = 0.0019 ± 0.0002 and m = 4.0 ± 0.2 for set #1; 

• Q, (C) AS m = 0.0022 ± 0.0002 and m = 3.8 ± 0.2 for set #2; 

• n (C) AS m = 0.0018 ± 0.0002 and m = 4.7 ± 0.2for set #3. 

A. Diffusion standard deviation analysis of the fluctuations 

We evaluate the stochastic equivalence made in Eq. (|25|) between the residual part of 
the discrete model, which makes use of the KL-decomposition, and the residual part of the 
continuous model, which makes use of a linear approximation, in two steps. Step 1 compares 
the size of the increments of the correspondent residual parts; and Step 2 adopts the standard 
deviation analysis (SDA) which is a statistical formalism to study the long-time correlation 
in a fractal time series. 

Because 0(t) = residual part, the increments are given by 9 t = 0(t) — Q(t — 1). We 
calculate the standard deviation, <jg of the increments {8 t } for each residual component 
estimated by means of the KL-decomposition and of the linear approximation respectively. 
Finally we calculate the average of the standard deviation, (erg), between the sixty <jg for 
each of the three cases. The results shown in Table I demonstrate the compatibility of the 
increments obtained with the residual parts of the KL-decomposition and the residual part 
of the continuous model. 

Now, let us suppose that a generic residual curve is given by the function O(t), see Eq. 
()25[) . that in this specific case is a kind of random walk around the ballistic part of the 
signal, which is the principal component of the KL-decomposition or the linear component 
of the continuous model. The SDA determines the scaling of the standard deviation of the 
diffusion process defined as 




N-T 



Q(t + r) - Q(t) - Q(t + r) - Q(t) 



(26) 



N-t-1 
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where 



N-T 



Q(t + r)-Q(t) 

N-T 



e(t + r)-e(0 = £ 



(27) 



i=0 



A" is the number of data points and the times t and r are measured in cycle period units. 
It is easy to realize that Eq. (|27|) ensures that £)(r = 1) = 1. In the presence of fractal 
statistics we would have, based on the discussion of anomalous diffusion in the Introduction, 



Figure 4 shows the SDA for the residual part of the KL-decomposition. Each set of graphs 
concerning the same crack data look quite similar. All three sets of graphs show that the 
curves have a initial scaling exponent approximately within the range 0.5 < (3 < 0.9. The 
mean curve value is represented by the curves with black circles in Figure 4. These early time 
values of /?, interpreted in terms of the random walks discussed in the Introduction, indicate 
that the residual parts of the signal manifest a persistent behavior, that is, a persistent 
correlation that lasts at least 10 consecutive cycles on average. 

For 10 < t < 100 the data presents a slight antipersistency with 0.4 < H < 0.5. Conse- 
quently, the residual process is initially strongly persistent, but asymptotically it is almost 
random. We observe that for 10 < r < 100 the mean scaling exponent is approximately 
H = 0.45 in the case of the linear continuous model and this is slightly larger than the 
scaling exponent in the KL discrete decomposition. This change in scaling is due to the fact 
that the principal part obtained with the KL decomposition extracts more information from 
the original signal than does the simple linear approximation. 

In the introduction we have explained that an initial anomalous diffusion that last up to a 
certain r as detected by Eq. could also be an artifact related not to some autocorrelation 
pattern in the data but to the transition from the initial geometrical properties of the 
distribution of the events {£,;} of a time series to the Gaussian shape of the asymptotic 
diffusion distribution. To check that the persistent behavior for r < 10 observed in the plots 
of Figure 4 expresses real correlation patterns, we repeat SDA of the data after randomizing 
the time series of the increments That is, for each crack data first we have the sequence 
{6t} defined as 9 t = 0(£) — Q(t — 1), then we shuffle {$t} and obtain a new sequence {0' t } 
and generate a new walk Q'(t) = Yll=i ^ an d finally we apply SDA to the new curve Q'(t). 
Figure 5 shows the result for the crack set #1 where the residual part is estimated with 
the KL-decomposition; for the other datasets the results are similar. Figure 5 clearly shows 



D (r) oc t p = r 



.a/2 



(28) 



14 



that after shuffling of the temporal order of the single increments {Ot}, the SDA of the new 
sequence gives a scaling value of approximately H = 0.5 and the persistent behavior for 
t < 10 observed in Figures 4 is absent. Thus, we conclude that the persistent behavior for 
r < 10 observed in Figures 4 expresses real correlation patterns in the fluctuations of crack 
growth. 

Figure 6 also shows that the the distributions of the scaling exponent seems to be quite 
uniform in the interval 0.5 < (3 < 0.9 (with a probability P > 0.9%) or, perhaps, as Figure 
6c shows better, there might be a slight prominence or skewness in favor of small value of 
(3. In any case, all figures show that the distribution of the scaling exponent for the residual 
components of the curve obtained with the KL-decomposition or the linear component of 
the continuous model practically coincide for all three datasets. This equivalence suggests 
that the continuous linear model essentially captures not only the dominant properties of 
the signal, as obtained through the KL-decomposition, see Eq. (J23j) . but also the stochastic 
properties of the residual signal, as suggested in Eq. (}25|) . 



B. Statistics of damage measure estimates 

We investigate the statistics of the damage measure estimates using a lognormal dis- 
tribution. This is in keeping with the analysis of several investigators who assumed the 
crack growth rate in ductile alloys is lognormal-distributed, see for example, the cita- 
tions in Sobczyk and Spencer jj. Other investigators have treated the crack length as 
being lognormal-distributed [stJ, rather than the residual fluctuations. The results of KL- 
decomposition in Eqs. (|16|) to (J19|) are consistent with these assumptions because Qo, which 
dominates the random behavior of fatigue crack growth, can be considered to be a perfectly 
correlated (ballistic) random process, whereas the non-negative, multiplicative uncertainty 
term (t, to) is a weakly (positively) correlated random process. Yang and Manning [39] 
have presented an empirical second-order approximation to crack growth by postulating a 
lognormal distribution of a parameter that does not bear any physical relationship to AS 
but is, to some extent, similar to Qq {AS) in the present model. 

Figure 7 shows the histogram of the slopes AS m Qq of the curves according to the 
continuous model for the experimental data presented by Eq. ()25j) . such as those shown in 
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Figure 3c. The histograms are fitted with the lognormal distribution p(x,/j,,a): 



1 



x\f2na 2 



cxp 



(ln(x) — ji) 
2^ 



2n 



(29) 



The measured parameters \i and a are recorded in Table II. Finally, the parameters /i and 
u are function of /iq = (x) and <7 2 = ((x — /iq) 2 ) as follows: 



/i = ln(// ) - cr 2 / 2 



and 



a 2 =ln 



(30) 



(31) 



Since the random parameter AS m Qq is not explicitly dependent on time, its expected value 
is obtained from Eq. (jl3|) as: 



t-t 



(32) 



which is readily determined from the ensemble average estimate from each of the data sets. 
Asymptotically in time we find for the variance of AS m Qo 

, 1 2\ 



a 2 = ((AS m fio-/xo) 2 ) 



ip (t, t 



t-t 



2 

- 



(33) 



so that the variance can be determined directly from the ensemble average estimate from 
each of the data sets. 



V. CRACK MODEL SIMULATION 

This section presents the results of Monte Carlo simulation of the fatigue crack damage 
process based on the model as it emerges from the stochastic analysis made in the previous 
section. The model that we introduce approximately reproduces the stochastic properties of 
both the ballistic or principal part of the fatigue crack growth and the associated fluctuations 
around it. The model consists in generating independently the fluctuation and the principal 
part of the fatigue crack damage in such a way they are statistically equivalent to the 
correspondent observations and then combining them. The crack model simulation is based 
on four steps: 



1(3 



Principal part or ballistic growth: We generate 60 values AS m Qo, lognormally dis- 
tributed according to Eq. where the parameters \i and a are given by the actual 
fit of the phenomenological distribution shown in Figure 6 and recorded in Table II. 
A sample of the curves AS m Qo(t — to) simulating the dataset #1 is shown in Figure 
8b. 

Residual part or fluctuations around the ballistic growth: We generate 60 fractal 
gaussian noise sequences {8[} each of length N of the original time sequence and with 
scaling exponent uniformly distributed in the interval 0.5 < (3 < 0.9. The standard 
deviation of each sequence is set equal to the mean standard deviation of the increments 
of the residual component of the data reported in Table I. To simulate the change of 
scaling exponent from persistent (for r < 10) to almost random (for r > 10), we 
section each fractal time series {8' t } into segments of length 10 within which the data 
would conserve the correlation, and finally we shuffle the position of these segments in 
the time series to reproduce a new time series {0 t }. These new time series will have 
persistent correlation for r < 10 and uncorrelated randomness for r > 10. Finally, the 
curve 0(t) is obtained by integrating the new sequence {8 t } and by detrending from 
it its linear component because the curve 0(t) is supposed having a zero mean. The 
SDA sample data analysis of an example of these synthetic residual data simulating 
the data set #1 is shown in Figure 8a. 

The ballistic growth estimated in the principal part and the associated fluctuations 
the residual part are combined according to Eq. (|THj) to obtain a simulated damage 
increment measure ip(t,to) for all sixty sequences and for the three datasets. Figure 
8c shows the simulated damage increment measure ip(t,t ) simulating the dataset #1. 

Finally, by using the respective value of the exponent m, reported in Section IIV[ for 
the dataset #1 and an one-dimensional root-finding computer algorithm, Eq. ([T2J) 
is inverted to obtain a simulated normalized crack length growth curves c t , as seen 
in Figure 8d. The similitude between Figure 8d and Figure la is noteworthy and 
the figures for the other data sets look qualitatively very similar; hence they are not 
presented in this paper. 



17 



VI. SUMMARY AND CONCLUSIONS 



This paper presents a stochastic measure of fatigue crack damage. We have focused on 
the correlation properties of the fluctuations around fatigue crack growth in ductile alloys. 
The model of crack damage measure indicates that the fluctuations around fatigue crack 
growth present strong correlation patterns within short time scales and are uncorrelated for 
larger time scales. These findings suggest that the random stochastic models adopted in 
the present literature for describing the crack growth dynamics should be augmented with 
short-time correlated stochastic models. 

The damage measure is modeled as an anomalous diffusion process that is obtained as 
a continuous function of the current crack length and of the initial crack length. Perhaps, 
the randomness in the damage measure estimate accrues primarily from manufacturing 
uncertainties such as defects generated during machining operations because such macro- 
defects are expected to drive the ballistic growth of cracks. This randomness is captured 
by a single lognormal-distributed random variable. Instead, the resulting diffusion process 
of estimated fluctuations around the ballistic growth of fatigue cracks is probably due to 
the inhomogeneity in the structural material because it is primarily associated with the 
micro-structure of the material, and is represented by a non-stationary fractional Brownian 
motion model. This non-stationarity manifests itself in the two scaling exponents occurring 
at different scales. Specifically, we observe a clear transition in the standard deviation 
analysis from an early time slope representing a strong persistence, (3 ~ 0.7 lasting for 
approximately r ?a 10 to a different slope asymptotically in time representing randomness, 
(3 ~ 0.5. This transition occurring at r ps 10 from a scaling regime to another indicates the 
scale at which a structure change of the ductile alloys occurs. 

The constitutive equation of the damage measure is based on the physics of fracture 
mechanics and is validated by KL-decomposition of fatigue test data for 7075-T6 aluminum 
alloys at different levels of (constant-amplitude) cyclic load. The damage estimate is 
statistically orthogonal to the resulting zero-mean estimation error in the Hilbert space 
L 2 (P) defined by the probability measure of the stochastic damage measure. As such, the 
damage estimate is often viewed as a best least-square linear estimate. However, we find 
that the KL-decomposition is statistically equivalent to the linear approximation in the 
continuum model that can be then used to simulate the fatigue crack growth in ductile alloys. 



18 



Acknowledgments 

The authors are grateful to Professor H. Ghonem of University of Rhode Island for provid- 
ing the test data of random fatigue crack growth. The work reported in this paper has been 
supported in part by the Army Research Office under Grant No. DAAD190110640. The 
first author thanks the Army Research Office for the support under grant DAAG5598D0002. 



Appendix: A Supporting Lemma 

Lemma: Let A (0 and B (£) be second-order real random variables; x((,t) and y(C,t) 
be zero-mean mean-square continuous (possibly non-separable) real random processes; and 
the real g (t) be almost everywhere continuous on an interval A such that, for alH e A, the 
following conditions hold: 

(i) A (C) = 5(C); 

(ii) (A (0 x (C, t)) = and (B (() y (C, t)) = 0. 
Then, the following mean-square identity 

A(Og(t) + x(C,t) = B(C)g(t) + y(C,t) 



yields 

x{C,t)=y{C,t) 
(A(()y((,t)) = \ Vt€ A 

(fl(0*(C,*)> = o, 

Proof: It follows from the above mean-square identity that 



Var [{A (0 -B(C)}g (t) + {x ((, t) - y (C, t)}] = 
which may be expanded to yield: 

Var [A (C) - B (()} g (tf + Var [x (C, t) - y (C, *)] 
+ ({A (0 - 5 (0} {x (C, t) - y (C, t)}) g(t)=0 

A combination of Condition (i) and Schwarz inequality yields: 



Var [x((,t)-y{(,t)}=0 
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and the remaining two identities follow from Condition (ii). 
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TABLE I: Values of the fiting parameters [i and a of the lognormal distribution (|29l) of the his- 
tograms shown in Figures 7. 





A* 


a 


Set #1 


0.58±0.05 


0.20±0.02 


Set #2 


0.74±0.05 


0.16±0.02 


Set #3 


0.42±0.05 


0.45±0.04 



TABLE II: Mean standard deviation of the increments of the residual part obtained with the K-L 
decomposition and the residual part of the continuous model. 





K-L 


Linear Model 


Set #1 : (a e ) = 


0.0024±0.001 


0.0025±0.001 


Set #2 : (a e ) = 


0.0024±0.001 


0.0025±0.001 


Set #3 : (a e ) = 


0.0038±0.003 


0.0043±0.003 
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FIG. 1: Experimental data of 7075-T6 aluminum alloy, [a] R=0.6 and Max stress=70.65 MPa; [b] 
R=0.6 and Max stress=69.00 MPa; [c] R=0.6 and Max stress=47.09 MPa. 
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[a] crack data 1 




1200 



1200 



200 400 600 800 1000 1200 
Time in units of 50 cycles (1 cycle = 100 ms) 



FIG. 2: [a] Curves ip(t) obtained with Eq. (|12|) for the experimental data of 7075-T6 aluminum 
alloy for set #1. The value of m used is m = 4.0. [b] Principal part of the K-L decomposition 
against [c] the linear approximation of the continuous model made in Eqs. (|23|) plus (|24|) of the 
curves tp(t). 
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FIG. 3: Increments Set against crack length q fit with Eq. (jHJ) (solid curve) for the experimental 
data of 7075-T6 aluminum alloy for the set #1: n AS m = 0.0019 ± 0.0002 and m = 4.0 ± 0.2. 
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FIG. 4: SDA for the residual part of the KL-decomposition (left figures) against SDA for the 
residual part of linear approximation (right figures) of the continuous model. Note the scaling 
transition at r ~ 10 from H ~ 0.7 to H ~ 0.5 in both cases for all datasets. 
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(shuffled crack data 1) K-L 



Q 




T 

FIG. 5: SDA for the residual part of the KL-decomposition after shuffling of the increments {Of}. 
Note the random scaling of H ~ 0.5. The data refer to set #1 and the comparison has to be made 
with Fig. 4 (crack data 1) KL. 
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FIG. 6: Histogram and probability density function of the scaling exponent (3 estimated by fitting 
the the range 1 < r <= 10 for each curve shown in the plates of Figure 4. Each plate compares 
the distributions obtained with the K-L decomposition and the continuous linear model for each 
of the three crack data sets. Figure 7d shows the histograms of all data. 
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FIG. 7: Histogram of the quantities AS m Qq of the continuous model of the experimental data 
presented by Eq. (|25|). The histograms are fitted with a lognormal distribution p{x; fj>, a) shown in 
Eq. (J2SJ). 
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[d] Synthetic set #1 
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FIG. 8: Synthetic data of crack length for Set #1. [a] SDA of the residual component; [b] ballistic 
growth; [c] damage increment ip(t)\ [d] crack length, compare with Figure la. 
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